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RIGOROUS DERIVATION OF ACTIVE PLATE MODELS FOR 
THIN SHEETS OF NEMATIC ELASTOMERS 

VIRGINIA AGOSTINIANI AND ANTONIO DESIMONE 


Abstract. In the context of finite elasticity, we propose plate models describing the spontaneous bending of 
nematic elastomer thin films due to variations along the thickness of the nematic order parameters. Reduced 
energy functionals are deduced from a three-dimensional description of the system using rigorous dimension- 
reduction techniques, based on the theory of P-convergence. The two-dimensional models are nonlinear plate 
theories in which deviations from a characteristic target curvature tensor cost elastic energy. Moreover, the 
stored energy functional cannot be minimised to zero, thus revealing the presence of residual stresses, as ob¬ 
served in numerical simulations. The following three nematic textures are considered: splay-bend and twisted 
orientation of the nematic director, and uniform director perpendicular to the mid-plane of the film, with 
variable degree of nematic order along the thickness. These three textures realise three very different struc¬ 
tural models: one with only one stable spontaneously bent configuration, a bistable one with two oppositely 
curved configurations of minimal energy, and a shell with zero stiffness to twisting. 


1. Introduction 

The interest in designing objects whose shape can be controlled at will through the application of external 
stimuli is fuelling a renewed interest in questions at the interface between elasticity and geometry. Which 
shapes are accessible to elastic sheets through the prescription of non-euclidean metrics that model states 
of pre-stress or pre-stretch induced by phase-transitions, plastic deformations, or growth [22]? Besides their 
fundamental mathematical interest [11, 24], these questions are very relevant in biology (e.g., in morphogenesis 
where shape emerges from growth and remodelling processes) and engineering (e.g., for motion-planning 
problems in soft robotics and, more generally, for the design of bio-inspired structures with programmable 
shapes). 

A general paradigm to generate bending deformations in thin films is to induce non-constant strains through 
the thickness^. These can in turn be triggered by the spontaneous strains associated with a phase transfor¬ 
mation. An example is provided by strips of nematic elastomers in which specific textures of the director 
have been imprinted in the material at fabrication. The process relies on pouring a nematic liquid between 
two plates which have been treated to induce a given uniform alignment of the director on one of them 
and a different one on the other one. This induces a non-constant director profile which is then frozen in 
the material by the photo-polymerization process that transforms a liquid crystal into a nematic elastomer. 
When the isotropic-to-nematic phase transformation takes place, the spontaneous deformations associated 
with it induce differential expansions along the film thickness, and hence curvature of its mid-surface. We 
refer the interested reader to, e.g., [30, 31, 34] for more details about the preparation of such materials, and 
to [2, 6, 10, 12, 13, 14, 15, 18] and the references quoted therein for further information on the mathematical 
modelling of their interesting behaviour. 

Two-dimensional models (plate models) for the bending deformation of thin films made of active material 
have already been proposed in the literature. They account for bending deformations through a curvature 
tensor (second fundamental form of the deformed mid-surface). The bending energy penalizes deviations of 
the curvature from a characteristic target curvature arising from the spontaneous strains triggered by a phase 
transition. Expressions for these bending energies are typically postulated on the basis of symmetry arguments, 
or deduced formally from an ansatz on the displacement fields (Kirchoff-Love assumption). By contrast, in 


^Another route, which exploits Gauss’ Theorema Egregium, is to induce curved configurations through nematic director 
textures generating (spontaneous) strains that are constant along the thickness but variable in the in-plane direction in such a 
way to be incompatible with having zero Gaussian curvature [3, 4, 5, 8, 25, 26, 27, 28, 29]. 
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our approach two-dimensional energy densities and target curvatures are deduced from 3D elasticity, i.e., from 
those geometric and material parameters that are available to the material scientists synthesizing the material 
and shaping it into a thin film. 

In this paper, employing rigorous dimension reduction techniques based on the theory of F-convergence, 
and following [32] , we derive new models for the bending behavior of thin films made of nematic elastomers in 
the regime of large deformations. Starting from three-dimensional finite elasticity, and considering the limit 
of a vanishingly small thickness, we obtain the following two-dimensional reduced “energy” functional 




(y) = ^ / { ^y{x) - A +j(^trAy(x')-trA^ | dx -b > 


( 1 . 1 ) 


whenever y is an isometry mapping uj (the planar domain representing the reference configuration of the mid¬ 
surface of the film) into Specific expressions of the 2D limit energy in terms of the parameters typically 
used in the theory of plates (such as the plate bending modulus) are given given in the right-hand-sides 
of (3.22) and (3.31). In (1.1) above, which is an expression of the type proposed in [3, 22, 33] to model 
the shaping of elastic sheets or of biological tissues, the symbol Ay denotes the curvature tensor, namely, 
the second fundamental form associated with the deformed configuration y{uj) and the coefhcients y, and 7 
are positive constants (material parameters characterising the three-dimensional stored energy density of the 
material). Moreover, the 2x2 symmetric matrix A is the target curvature tensor and e is a nonnegative 
constant. The characteristic quantities A and e are deduced from the three-dimensional model and given 
by explicit formulas, issuing from the specific variation of the spontaneous strain along the thickness. The 
constant e is irrelevant in the selection of energy minimising or equilibrium shapes y. However, a term e > 0 
is typical for those cases in which the spontaneous strains of the 3D model are not kinematically compatible^. 
Thus, just like its parent 3D energy functional, the limit 2D energy (1.1) can never be minimised to zero: 
there will always be energy trapped in the system, indicating the presence of residual stresses. 

Two special geometries of the director field are of particular interest, since they have been realized in 
practice in the laboratory. In the splay bend geometry (SB), also called hybrid in [30, 34], the explicit 
formulas we obtain for A and e are 

= ^Sb(<5o) = ^^diag(- 1 , 0 ) and e = esB{5o) = + ^)5l ^ [a;], 


where 5q is a positive constant, with dimension of inverse length, which quantifies the variability along the 
thickness of the spontaneous strain (see (2.3)-(2.4) and the first formula in (3.15)). We recall that, in a 
three-dimensional film with splay-bend geometry, the director continuously rotates by 7r/2 from planar to 
vertical alignment (see (2.7) and Figure 1). The other geometry we consider is the twisted one (T), see [31] 
and [34], where instead the director (continuously) rotates perpendicularly to the vertical axis from a typical 
orientation at the bottom of the film to another typical orientation at the top of the film (see (2.7) and Figure 
1). In this case, it turns out that 


12do 


A = At(So) = diag( - 1, l) 


and 


e = eT(do) = 


- 47r2 - 48 




The difference in the formulas for the two cases arises because of a different distribution of spontaneous strains 
along the thickness, see (2.3)-(2.4) and (2.7). 

The two geometries of the director field described above lead to plates with two very different structural 
behaviours. Both of them arise from kinematically incompatible spontaneous strains, which generate residual 
stresses leading to a strictly positive constant e. They differ in the fact that in the splay-bend geometry, the 
integral term in ( 1 . 1 ) can be minimised to zero (by any developable surface y(uj) whose second fundamental 
form Ay coincides with Asb)- Hence, the target curvature Asb is the curvature the plate spontaneously 
exhibits in the absence of external loads (spontaneous curvature). By contrast, in the twisted case, also the 


^To put the incompatible nematic elastomer cases in perspective, we also analyze a kinematically compatible case where the 
three-dimensional spontaneous strain distribution along the thickness depends quadratically on the thickness variable. We show 
that, as expected, this distribution leads to plates with no residual stresses and where the limiting energy corresponding to (1.1) 
attains its minimum value zero. 
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minimum of the integral term is strictly positive. In fact, there exists no isometry y such that Ay = At, 
because in a developable surface the product of the principal curvatures must be zero at each point of the 
surface. This means that, in fact, the target curvature At is never observed in the absence of external loads. 

The spontaneous curvature exhibited in the absence of external loads by a nematic film with twisted texture 
cannot be read off directly from the target curvature, but is has to be computed by minimising the integrand 
in (1.1), subject to the isometry constraint. It turns out that this system has two distinct configurations of 
minimal energy, with opposite curvature, hence it is bistable. By contrast, in the splay-bend case, there is only 
one stable bent conhguration. Motivated by these observations, we also consider the following geometry for the 
nematic parameters: a uniform director orientation perpendicular to the mid-plane of the film, with variable 
degree of nematic order along the thickness. Even though this configuration has not yet been realised in the 
laboratory, it leads to a very interesting mechanical behaviour. Namely, a structure possessing a continuum 
of spontaneously bent, minimal energy configurations, representing a shell with zero stiffness to twisting. 

The rest of the paper is organised as follows. In Section 2, the 3D elasticity models are presented and 
a discussion of the kinematic compatibility of the 3D spontaneous strains is provided. Then, in Section 3, 
we present the theoretical basis for our dimension reduction procedure, and the derivation of the formulas 
allowing to deduce the target curvature A and the constant e from 3D elasticity. This is the content of 
Theorems 3.5, 3.7 and of formulas (3.22) and (3.31). As already mentioned, we work in the framework of 
the dimension reduction approach which traces back to the seminal paper [17]. In particular, to obtain our 
results we use the plate theory for stressed heterogeneous multilayers developed by Schmidt in [32], which 
has recently motivated some new computational schemes [7]. Our models are valid for arbitrarily large elastic 
deformations. A plate model covering the regime of small deformations has been presented in [21]. 

Section 4 is devoted to the physical interpretation of our results: We derive explicit formulas for the 
deformations realising the minimal free-energy of the (reduced) plate models, which represent the configuration 
the nematic sheets exhibit in the absence of applied loads. We show that there is one spontaneously bent 
configuration in the splay bend case, while there are two distinct ones in the twist case. Thus, twist nematic 
plates are bistable structures, a fact that has gone unnoticed until now, and has not yet been observed in 
the laboratory. Moreover, the behaviour of splay-bend and twisted nematic elastomer sheets is compared to 
the case in which the nematic director field is constant (perpendicular to the mid-surface), and the thickness- 
dependence of the spontaneous strain is induced by the variation of the degree of nematic order along the 
thickness. Although a system like this has not yet been synthesised in the laboratory, we hope that the 
predictions of our model will motivate researchers to investigate experimentally the mechanical response it 
would produce. In fact, our prediction is that, in the thin him limit, this texture should produce a plate with 
soft response to twisting, see Figure 3. 

2. SpLAY-BEND and twisted NEMATIC ELASTOMERS THIN SHEETS 

In this section, we present a three-dimensional model for a thin sheet of nematic elastomer with splay-bend 
and twisted distribution of the director along the thickness. The kinematic compatibility of the corresponding 
spontaneous strains is discussed in Subsection 2.2, where the case of strains distributed quadratically along 
the thickness is analyzed as well. 

2.1. A three-dimensional model. We consider a thin sheet of nematic elastomer occupying the reference 
conhguration 

nh=uj X {-h/2,h/2), (2.1) 

for some h > 0 small, where w is a bounded Lipschitz domain of with sufficiently regular boundary. 

Notation 2.1. Throughout the paper we will denote by { 01 , 62 , 63 } the canonical basis of and by z = 
(zi,Z 2 ,Z 3 ) an arbitrary point in the physical reference conhguration D/i. The term “physical” here and 
throughout the paper is used in contrast to the corresponding rescaled quantities we will introduce later on. 
Also, SO(3) is the set of the 3x3 rotations and I S SO(3) the identity matrix, whereas the symbol I 2 denotes 
the identity matrix of 

We suppose the sheet to be heterogeneous along the thickness with associated stored energy density 

: (-h/2, h/2) X M3x3 ^ [0, -hoo]. 
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More precisely, in the two models we are going to consider, the 23 -dependence of the energy density is induced 
via the 23 -dependence of the spontaneous strain distribution. 

If n S is a unit vector representing the local order of the nematic director, the (local) response of the 
nematic elastomer is encoded by a volume preserving spontaneous strain (technically, a right Cauchy-Green 
strain tensor) given by 

L{n) = a^wSin + a~3 (I — wSi n), (2.2) 

for some material parameter a > 1, which is usually temperature-dependent. Suppose that the nematic 
director n varies along the thickness according to a given function 23 1 —>■ n^{z 3 ) and coincides with two given 
constant directions at the top and at the bottom of the sheet: 

n^(—h/2) = rib, n^{h/2) = rit, for every small h > 0, 

for fixed nb,ut S S^. The through-the-thickness variation of the nematic director translates into a variation 
of the corresponding spontaneous strain according to ( 2 . 2 ), namely, 

Chizs) ■■= L{n^{z3)) = a^^^n''( 23 ) (g) 71 ^( 23 )-I-- n'‘( 23 ) (g) n'‘( 23 )). (2.3) 


Notice that, in this expression, we allow the material parameter a to be Ii-dependent. 
now on we will assume that 

1 ^ 

O/t — 1 + 

ho 


More precisely, from 


(2.4) 


where ag is a positive dimensionless parameter, while ho and h have the physical dimension of length. This 
assumption is easily understandable if one thinks that curvature is related to the ratio between the magnitude 
of the strain difference along the thickness and the thickness itself. Hence, the linear scaling in h in (2.4) is 
needed in order to obtain finite curvature in the limit /i —> 0. Observe that Ch{zo) is positive definite for every 
23 S (—/i/2,/i/2) and every h > 0 sufficiently small. 

In the framework of finite elasticity, a prototypical energy density : {—h/2, h/2)xM.^^^ —>■ [0, 00 ] mod¬ 
elling a nematic elastomer is 


w\zo,F) 


h 
2 L 
- 1-00 


iF^F)-cj;\zo)-3 


2 log(det F) 


W^oiidetF) 


if det F > 0 , 
if det F < 0 , 


(2.5) 


where /r > 0 is a material constant (shear modulus) and the function Wyoi '■ ( 0 , 00 ) —>• [ 0 ,oo) is around 1 
and fulfills the conditions: 


W.oiit) = 0 ^ t = I, W,oi{t) ^ 00 as t ^ 0+, 1T;„,(I) > 0. 

It is easy to show (see Remark 3.6, (3.16), and (3.24)) that is indeed nonnegative and such that 

u;^ 23 ,F) = 0 iff Fe SO(3 )Vc^(23). 

Expression (2.5) is a natural generalization, see [1], of the classical trace formula for nematic elastomers derived 
by Bladon, Terentjev and Warner [9], in the spirit of Flory’s work on polymer elasticity [16]. The presence 
of the purely volumetric term Wyoi (det F) guarantees that the Taylor expansion at order two of the density 
results in isotropic elasticity with two independent elastic constants (shear modulus and bulk modulus). 

If {fh}h> 0 : with represents a family of applied loads, the (physical) stored elastic energy 

and total energy of the system associated with a deformation u : —>■ are given by 

i\v) = I w^{zoyv{z))dz, = i^{v)-j h-vdz, ( 2 . 6 ) 

respectively. 

Let us now focus on the nematic director field in the splay-bend and twisted cases, which we denote by 
inlgg and n^, respectively. We recall that these distributions are solutions to the problem 

min / |Vnpd 2 , 

n{—h/2) = J 
n(h/2) = n* 
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where in the splay-bend case = ei and n* =63, whereas in the twisted case n*’ = ei and n* = 62- We have 
/eos(f + ff)\ /cos(f + ff)\ 

'^SbI^s) — I 0 1 ) B.rp^z^) — I sin(^ -I- f x) 1 > ^3 € (~^/2,/ 1 / 2 ), (2-7) 

vsin(f + ff)y V 0 j 

and we refer the reader to Figure 1 for a sketch of these two geometries. 

We define the (physical) spontaneous strain distributions Ch,SB and Ch,T as that in (2.3) with nJgg and 
in place of n^, respectively. Correspondingly, we denote by Wgg and the stored energy densities, by S’gg 
and S'^ the stored energy functional, and by and the total energies. 



Figure 1. 


Sketch of the splay-bend director field (left) and of the twisted director field (right). 


2.2. Kinematic compatibility. Here, we want to discuss the kinematic compatibility of some given field of 
(physical) spontaneous strains. Let O be the (physical) reference configuration of a given system and suppose 
that it is a simply connected open subset of We say that a smooth map G : O —> representing a 

distribution of spontaneous strains and such that G{z) € Psym(3) for every z € O, is kinematically compatible 
if there exists a smooth function v : O —representing a deformation and such that det Vn( 2 ;) > 0 for 
every z € O, satisfying 

= G in O. (2.8) 

Following [11], we reformulate this concept in the framework of Riemannian geometry. In order to do this, 
let us denote by 5^3 the Euclidean metric of and recall that for a given immersion ip : O —(K^, ) the 
pull-back metric of via p is the metric h defined in O by the identity 

h\p{X,Y) = (d(/ 5 |p[X],d(p|p[F]) , for every X,Y e TpO = 

The pull-back metric of g via p is usually denoted by p* g\ip(o)- If (•z*)i=i, 2,3 and (a:“)Q=i^ 2,3 are systems of 
coordinates for O and p^O), respectively, the above identity specialized to X = d/dz]^^ and Y = djdz^^^ gives 

_ jj3 dp°^ dp^ 

az* Ip Ip’ 

where p'^ := a;“ op. If, in addition, we assume (a;“)ct=i, 2,3 to be the standard Euclidean coordinates, then the 
coefficient is just Note that here and in what follows the Einstein summation convention for the 

sum over repeated indices is adopted. We identify the given spontaneous strain distribution G with a metric 
defined in O, so that asking if there is an orientation-preserving deformation v : O —>■ such that ( 2 . 8 ) 

holds true corresponds to seeking for a local diffeomorphism v : O —such that the pull-back metric of 
gs.3 via v coincides with the metric G. In formulas, 

= "a/3 


dz^ Ip dzi |p’ 
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where we have fixed standard Euclidean coordinates in the target manifold v{0). Note that since w is a local 
differmorphism, the equivalence 

^*550)=^ in (2-9) 

establishes a local isometry (through v) between the manifolds {0,G) and Now, we have 

from Theorem 1.5-1 and Theorem 1.6-1 in [11] that, since O is simply connected, a necessary and sufficient 
condition for (2.9) to hold is that 

Riemc = 0 in O, ( 2 - 10 ) 

where Riemg is the fourth-order Riemann curvature tensor associated with the metric G. We recall that, in 
the given local chart (- 2 :*)i=i, 2,3 of O, the (3, l)-coefficients R^^^’s of Riemc = Rb^(d 2 :* 0 dz^ 0 dz^ 0 djdz^') 
are given by 

3 d 

Til __ ^ -pi _ -pi I -pi -pS _ -pi pS 

Wj/C ■“ Qz 3 ij ^ ik ^ks^ijt 

where the Christoffel’s symbols T^^ ’s are defined as 

:= r,,, := i (a,G,7 + O^Gu - diG^j ), (2.11) 

and the symbols G^-^ ’s stand for the components of the inverse G~^ of G. To simplify the computations it is 
sometimes useful to introduce the (4, 0)-coefficients RiijkS of Riem^, defined as 

Riijk ■= = djViki — dk^iji+T^-Ti;ip+T^iJ'jip. 

It is clear that Riemg = 0 if and only if Rajk = 0 for every I, i,j, k G {1, 2, 3}. Finally, let us recall that, since 
we are in dimension 3, condition (2.10) is equivalent to Ricc = 0 in G, where Ricq denotes the second-order 
Ricci curvature tensor associated with G, which is defined as Rice = Rij dz^ 0 dz^ , with 

R., := dir‘^ - djT‘, + 

From now on in this section, we restrict our attention to the case where O = flh (see (2.1)) and the spontaneous 
strain distribution G is a function of the thickness variable Z 3 G (—/i/2, h/2). Note that a material point of fl/j, 
normally referred to as a point of components (zi, Z 2 , Z 3 ) throughout the paper, is a point of the manifold Qfi 
with coordinates (z^, z^, z^) from the point of view of Riemannian geometry. In the following subsections, we 
discuss the kinematic compatibility of Z 3 i-A- G(z 3 ) in three cases: the case where G(z 3 ) depends quadratically 
on Z 3 and two cases (splay-bend and twisted nematic elastomer sheets) where the dependence of G{z 3 ) on Z 3 
is more complicated and gives rise to incompatible strains. Throughout this section we use the variable t in 
place of Z 3 and we use the index/apex “t” in place of “3”. 


t G {-hj2,hj2), 


( 2 . 12 ) 


2.2.1. The splay-bend case. In this case, setting 

^ / X 'TT tt 

m ■-= j + 

and looking at (2.3) and (2.7), we have that, up to a multiplicative constant, the spontaneous strain distri¬ 
bution is given by 

1-f (a^-l)cos2//. 0 

0 1 

sin(2A) 0 


G — G{t) — I + {ah — Gj'nJgB ® — 


(^)sin(2A) 

0 

l + {ah- 1) cos^ fh 


G-^(t) =1 + 


1 


- 1 




’ ’^SB — 


( 0 (^)sin(2A) \ 

0 10 

V (^)sin(2A) 0 l+(^-l)cos2//. ^ 


:= d,r\,-d,T\,+T\,T^,,+Tlr^B = - 5 *(r]i+r 2 ,)+rl,rJ,+r[,rJ, = (2.13) 


It turns out that the coefficient Rit of Ricc has the quite simple expression 

Rlt ■= O’ii It-Dti liWJ- IJ-TJ. (fcl- i; 

where we have first used the fact that the Christoffel symbols depend only on t and secondly the property 

whenever 2 G {i,j,k}. 


r'^- = 0 
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This can be easily checked using the definition of T^^ in (2.11). The same definition and simple computations 
also give 


ill — 


^ It — 


pi — 
-L ++ — 


—— — 


(a — 1)^ / 77 


(^) ®“VliC0sVli, 


-(a - 1) sin fh cos fh 
{a - 1) (cos^ fh - sin^ fh) 


1 + ( - - 1 ) cos fh 


1 + ( - - 1) sin fh 


(a-1)^ 


(^) fhCOs"^ fh, 


pt _ 
i 11 — 




(a - 1) sin fh cos fh 


1 + ( - - 1) sin^ fh 


Plugging these expressions into (2.13) yields 

Pit = sin fh cos fhicos'^ fh - sin^ fh) = sin . 

Thus, we can conclude that is not identically zero in In turn, Rico is not identically zero, so that the 
splay-bend spontaneous strain distribution is not kinematically compatible. 


2.2.2. The twisted case. In this case, following the same notation as in (2.12), we have 

/ 1 -b (a/j - 1 ) cos^ fh sin{2fh) 0 \ 

G = G(t) = I-b (a/^ - l)ny (g) = I (^^) sin(2/,i) 1 + {ah - 1) cos'^ fh 0 , 

0 0 1 / 

(^)sin(2A) 0^ 

(^)sin(2A) l+(^-l)cosV/. 0 

0 0 1/ 

For the twisted geometry, the coefficient Rtt of Ricc has a simple expression. Indeed, we have 


G ^ (t) = I -b-1 (g) = 


a-h 




R„ := diTi, - atT: 


■p/ Y\k 

Ik^ tt ' 


■pZ p'/c 

^ tk^ tl 


= -9t(r( 




■^tA) = 




+ 2Tl,Tl, + {Tl,)^ 
(2.14) 


since 

rtt = 0 for every fc = l,2,f, and TJ^ = -T^^ = - sincos//,. (2.15) 

This can be easily checked from the definition of the Christoffel symbols in (2.11). Similar computations yield 



Using these formulas together with (2.14) and the second equation in (2.15) gives 


Pit 


{a - 1)^ 

2a \2h) 


This fact implies, in particular, that Picc is not identically zero and in turn that the twisted spontaneous 
strain distribution is not kinematically compatible. 


2.2.3. The quadratic case. In this subsection, we consider the case where 

G = G{t) = l + tA + i^ B, 


t G (-/i/2, V2) 


(2.16) 
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for some diagonal matrices A = diag(Aii, ^422, Att) and B — diag(Bii, B 22 , Bu)- Note that G{t) =: diag(Gii(<), G 22 {t),Gtt{t)] 
is positive definite for every t sufficiently small. Elementary computations yield 



Rii — — 

R 22 = — 

^ “4 1 ^dt \^Gii 
Rij =0, for every i ^ j, 

where Ga is the derivative of Gu with respect to t. Now, set 

^:=logGii, ?7:=logG22, 

so that Gii/Gii = G22/G22 = f], GttIGtt = r, and in turn 



r logGtt, 






+ r] — T 


R22 - -4 ( - 


(^) [^V + v{i + V-r) 


Rtt — ~ 


1 

4 L 


2{i+ii)-f{i + fi) + Uf + {fif 


The condition Rice = 0, which guarantees the kinematic compatibility of t 1 —>■ G{t) as discussed above, is 
then equivalent to the following system of ODEs: 

f 21 + ^ + 77 - f j = 0, 
f 277 + 77^^ + 77-7-^ =0, 

i 2{^ + ij)-T{i + ^) + {ir + {f,r=0. 


Solving this system translates into compatibility conditions on A and B in ( 2 . 16 ). It then turns out that a 
spontaneous strain distribution t 1—>■ G{t) of the form ( 2 . 16 ) is kinematically compatible if and only if one of 
the following four conditions is satisfied: 

(*) ^11 = 4I22 = Att = 0 and Bn = B22 = Btt = 0 , 

{a) A22 = Att = 0, B22 = Btt = 0 , and Bn = An/d ^ 0 , 

(m) An = Att = 0 , Bn = Bu = 0 , and B22 = ^22/4 ^ 0 , 

(iv) An = A22 = 0 , Bn = B22 = 0 , and A^t + Bft ^ 0 . 

Note that the first condition corresponds to the trivial case G = I and the second one tells us in particular 
that a strain of the form 

/ l — 2kt + k^t‘^ 0 0\ 

G(t) := 0 10 , t€{-hl2,hl2), ( 2 . 17 ) 

V 0 0 1/ 

for some constant fc S M \ { 0 }, is kinematically compatible. A prototypical deformation v giving rise to such 
G can be provided in the following way. Let I be an open interval of M, let 7 : / — y be a smooth curve, 
and define 

r(s):=7'(s), Ar(s):=^^, B{s) := T{s) A N{s\ 
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for every s G /, where the apex stands for differentiation with respect to s. Suppose that the curve is 
parameterized by arc length, so that \T\ = 1 and the curvature k is defied as k := \T'\. Then the Frenet- 
Serret formulas read 

{ T' = kN, 

N' = -kT + rB, 

B' = -tN. 


Note that multiplying the first equation by N gives k = T' ■ N = —N' ■ T. Let us restrict to the case of B 
being constantly equal to 62 , where { 61 , 62 , 62 } is the canonical basis of This means that 7 is a planar 
curve and the above formulas imply in particular that r = —B' ■ N = 0 and = k'^. Now, let us define 

V : flh —> as 

v{s,Z2,t) := 7(s)+tfV(s)+ Z2e2, 

where we have supposed ujh = ui x (—/i/2,/i/2) with ui = I x J, for some open intervals /, J C R. Then 
Vv = {T + t N' \ e 2 \N) and therefore 

/ |r|2 + 2 tT-iV' + t2|iV'|2 e 2 -(T + tiV') N\T + tN')\ / l-2kt + k^f 0 0\ 

Vu^Vu= e2iT + tN') |e2|2 e2 ■ N = 0 10. 

\ N-lr + tN') e2-N |7V|2 J \ 0 0 1/ 

Supposing the curvature k to be constant, we have thus derived a strain G of the form (2.17). 

Finally, note that the analysis performed in this session shows that in the linear case where G{t) is of the 
form 

G = G{t)=l + tA, te{-hl 2 ,hl 2 ), 

for some diagonal matrix A ^ Q, the kinematic compatibility of the spontaneous strain distribution is never 
fulfilled. 


3. Derivation of the plate model 

In this section, we first rewrite the three-dimensional model previously introduced in a rescaled reference 
configuration. Then, in Subsection 3.2, we recall two rigorous dimension reduction results of compactness 
and F-convergence. This mathematical technique is subsequently employed in Subsection 3.3, where our main 
results. Theorems 3.5 and 3.7, are stated and proved. 


3.1. The rescaled three-dimensional model. As it is standard for dimension reduction techniques, let us 
now operate a change of variables in order to rewrite the energies in a fixed, h-independent rescaled reference 
configuration. 

Notation 3.1. We denote by a; = (a;i,a; 2 ,X 3 ) = {x' ,Xz) an arbitrary point in the rescaled reference conhgu- 
ration fl := uj x (— 1 / 2 , 1 / 2 ). 

For every h> 0 small, we define the rescaled energy density : (—1/2,1/2) x R^^^ —[0, -l-oo] and the 
rescaled applied loads fh ■ R^ as 

W^{x 3 ,F) := w'"(hx 3 ,F), fh{x) := fh{x ^hx^). (3.1) 

Note that fulfills 


W’^[x3,F)=Q 


Setting 


^hV ■= 


dx2y 


iff 

dx^y 


F &SO{i)JGh{z3), 




dx^y 

h 


Ghixz) ■■= Chihxs). 


for every y : Q 


(3.2) 


the correspondence between the original quantities and the rescaled ones is through the formulas 
i^{v) = v{z) := y{z', z^/h) a.e. z in flh. 


(3.3) 
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Here, the rescaled stored elastic energy functional and the rescaled total energy functional are defined, 
on a deformation y : fl —>■ K^, as 


J W^ix3,Whyix))dx, ■^^{y) ■= <^^iy) - j h ■ ydx. 


(3.4) 


n n 

Following the notation already introduced in Section 2.1, we use the indexes SB and T to denote the quantities 
related to the splay-bend case and twisted case, respectively. Hence, we write Ch,SB, ’^sb 

for the splay-bend model, and Ch,T, and for the twisted model. 

We now focus attention on the (rescaled) spontaneous strains Chix^). Looking at (2.7), we first note that 
for both models n^{hx 3 ) is independent of h, namely 


cos I 


NsBixs) := ngsihx^) = 


+ ^Xs) 


0 

-I- 

.4 ^ 2 ‘ 


sm (-I- fxa) 


Nrixs) := n^^kxs) = 


cos (f -h f X 3 ) 

sin(f + 

0 


(3.5) 


for every 313 G (—1/2,1/2). Hence, referring to the (above) definition of Ch and to expression (2.3), we have 
for the splay-bend case as well as for the twisted case 

Chixs) = af/^N{x3) (g) N{x3) + - N{x3) g) N{x3)) 


f 2/3 -l/3\ 

= [o-h - o-h ) 


a/i - 1 


■ N{x3) g N{x3) 


aoh 


= I + w— N{x3) g N{x3) - - + R {X 3 ), 


where ||i?^||oo = o{h) and || • ||oo is the norm in the space L°°((—1/2,1/2),IR^^^) 
equality we have plugged in expression (2.4) for au and used the expansion 

2 


(3.6) 

Note that in the third 


2/3 - 1/3 Oi^h 1 / a^h 

1. — n. ' = -— - 


hn 


3 \ ho 


o{h? 


3.2. A rigorous mathematical result for the limiting theory. For the convenience of the reader, we 
collect in this section, in a slightly simplified version, two results proved in [32] (Theorems 3.3 and 3.4 
below), which we are going to use later on. In this paper, an arbitrary family of energy densities : 
(—1/2,1/2) X [0, -l-oo] is considered, with the property that 

W'^ixo.F) = Wo{F{l + hB^{x3))). (3.7) 

where the function Wq : [0, -l-oo] satisfies Assumption 3.2 below, and 

B^ ^ B in L“((-1/2,1/2),K3 x3 )^ 10. 

For each small h, let us introduce the functional : W^’^(n,K^) —>■ [0, -l-oo], defined as 

<^^2/) := I W’^{x3,^hy)dx, 
n 

with V/i given by (3.2). Recall that here and throughout the paper = uj x (—1/2,1/2) and w C is a 
bounded Lipschitz domain with sufficiently regular boundary. More precisely, for the following theorems to 
hold, it is required that there exists a closed subset E C 9a; with = 0 such that the outer unit normal 

exists and is continuous on duj \ E. 

Assumption 3.2. The function Wq : —>■ [0, -l-oo] fulfills the following conditions: 

(i) it is in a neighborhood o/SO(3), and it is minimised at I; 

(ii) it is frame-indifferent, i.e. Wo{F) = Wo{RF) for every R G SO(3). 

(iii) there exists a constant C > 0 such that for every F G 

Wo(F) > Cdist^(F,SO(3)). 
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The following result states that a sequence {yh} which bounds the energy by a factor h? converges 
(up to subsequences) to a limit that is constrained to the class of isometric immersions of uj into the 

three-dimensional Euclidean space, namely 

Aso := [y G W2’2(a;,K3) : {VyfVy = h a.e. in w}. (3.8) 

Theorem 3.3 (Compactness). If {y^} C W^’^(n,R^) is a sequence such that 

I W>^{x3,Vhyh)dx<Ch‘^ (3.9) 

n 

for every h> 0 small, then there exists a (not relabelled) subsequence such that 

Vhy^ — >{Vy\v), as hfO, strongly in L^(n,R^^^). 

Moreover, the function x i—>■ (V'y | u) belongs to W^’^(r2, is independent of Xs, and 

[y'y I v) {x') G SO(3) for a.e. x' G uj. 


Before proceeding, let us introduce some more notation and denote by M G R^^^, the quadratic 

form D^bFo(I)[M, M], where D^Wo(I) stands for the second differential of Wq evaluated at I. Moreover, define, 
for every G G R^^^, 


02 (G) := 


min O 3 


G 

b 

0 

a 


(3.10) 


and in turn 


1/2 

02(G) := J Q2{D + tG + B{t))dt, 

-1/2 


(3.11) 


where B is obtained from B by omitting the last row and the last column. 


Theorem 3.4 (T-convergence). The functionals fh? T-converge as h fQ, with respect to the strong and 
the weak topology 0 /W^’^(n,R^), to 

^ J Q2{Ay{x'))dx' if y & Aiso, 

UJ 

-foo otherwise in W^’^(r2,R^), 

where Ay denotes the second fundamental form associated with the surface y(uj). 

Recall that the second fundamental form of y{uj) at a point y{x') can be expressed as (V'?/(x'))^W(a:'), 
where v := d^^y A d^^y- 




^'*'"( 2 /) := 


3.3. Splay-bend and t-wisted nematic elastomer plates. We want to apply the theory presented in the 
previous section to our two models. We first focus on the splay-bend case, whose associated rescaled stored 
energy density, considering expression (2.5) together with(3.1) and (3.6), is given, for every 0:3 S (—1/2,1/2) 
and every F G R^^^ with det F > 0, by 


WL(^3,F) 


T 

2 


{F^F) ■ C^^sB 


(xs) — 3 — 2 log(det F) 


+ lEi,o;(det F). 


Recall that 


Ch.ssixs) — I + 


ho 


(^Nsb{x3) C) Nssixs) 


+ H-SBi^s), 


WR’^sbIU = o{h) (3.12) 


(see (3.5) for the definition of Nsb)- Defining 


Wo{F) := 



-3-21og(detF) 


-h Wyoi{detF), 


(3.13) 
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for every F G 


3x3 


with det F > 0, and setting 


Uh,SB ■= y Ch,SB, 

yields Wgg{x 3 , F) = Wo[FU f^^gg{x 3 )). Note that 


(3.14) 


u-^ 


h,SB 


= 1 + h 


-Jo Msb - - + ^ 


r,h 

^SB 


with llfg^lloo = o{h), where we have used the notation 


Ci^n 

Jo := A7~) Msb ■= Nsb ® Nsb = 

ZtlQ 


cos^/i 0 |sin(2/i) 


0 


0 


0 


isin(2/i) 0 sinVi 


Here, the function 0:3 1—>■ fiixs) is defined as in (2.12), with h = 1 . All in all, we can write 


W^si^,,F) = Wo{F{l + hB^gs{x,))), 


Tjh ._ 

^SB ■ — 


—Jo ( Msb ~ a ) + 


r.h 

^SB 


(3.15) 


(3.16) 


Since Hr^sHoo = o{h), we have that —s> Bsb inL°“((-l/2,1/2), , where Hss := -Jq {Msb{x 3 ) - |). 

In turn, also in view of Remark 3.6 below, we have shown that the splay-bend model introduced in Section 2.1 
perfectly fits the mathematical theory summarized in the previous section. Hence, we have to compute the 2D 
energy density according to formula (3.11). First of all, we have that Q^IM) = 2fi |symMp -f tr^M. 

Using this expression, we can compute Q 2 for every G G (see (3.10)): 

+ lU";(l)(trG + a)4 


Q 2 {G) = min < 2/1 

aGR ^ 


symG 

b /2 

~ 1 FJT~ 

a 


= 2/i|symGp-k W'"„;(l)tr2G-kmin (2/i-k lU"o;(l))a^-k 2 VlU/;(l)tr G a = 2/i (|symGp-k ytr^G) , 

aGM L J 

(3.17) 

(3.18) 


having introduced the notation 


7 := 




2/i + lU",(l)- 


Finally (cfr (3.11)), note that Bsb is given by 
I 2 


Bsb = -Jo ( Msb “ y ) > 


with MsBixs) := (NsBixs) iSi Nsb{x 3 )) = ( 4 ^ 2 ^ 3 ) 0 


We are now in the position to compute, for every G G 
1/2 

Q2 Sb(G):= min [ Q 2 {D+ tG + ISsBit)) dt 


-1/2 


1/2, 


= 2/i min 

DGSym(2) 


- 1/2 


r 

D + tsymG- JoMsB(i) + yl2 


-I- ytr^ ( D -I- t G - Jo MsB^t) + ^l2 )/■ dt. 


The integrals 


1/2 

1/2 

1/2 

J \MsB\‘^dt = 

/ '“■‘(l + 5*)'^‘=5' 

/ 

-1/2 

-1/2 

-1/2 


1/2 

ttyMsBdt= J teas' 

-1/2 


(4 + 2/ 


dt = —- 


and other elementary computations imply that Q 2 g^{G)/{ 2 fi) equals 
^(isymGp+ytr^G) +^(symG-diag(l,0)+7trG) + 


<^0+ min USB (F), 
I)eSym(2) 
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where 


qsB{D) := \D\'^ +jtr'^D - (5o 


symD ■ diag(l,0) + 


D 


It is easy to see that 


min 

DeSym(2) 


qsB{D) = gsB(diag(do/6,-(5o/3)) = 


and in turn that 

Q2,sb{^) = 2 ^ 


1 

12 


(isymGp+ytr^c) + ^(symG-diag(l, 0 )+ 7 trG) + 


1 + 71 


It is again a simple computation showing that there exist constants asB, Psb G R and Asb G Sym(2) such 
that 

Q 2 ,SBiG) = asB Q 2 [G - Asb^ + PsB, for every G G Sym(2), (3.19) 

and they are given by 


OiSB = 


1 TT 12 do 

Asb = ^^diag(-l, 0 ), 


PsB = m (1 + 7)'5o( 


TT^ - 12 


(3.20) 


To state our result, let us define the functional ^ 5 ™ : A, 
(3.8), as 


—[0, 00), where Aiso is the class defined in 

(y) ■■=11 Q2.sBiM^'))d^' 

OJ 

= Ay(a;')-^diag(-1,0) +7 (h^(x') + ^) jda;'+ (1 + 7)^0 g (3-21) 

LO 

Here, the symbol Hy denotes the mean curvature of y{uj), hence Hy = trHy. Note that for every y S Aiso we 
have that |Ay| G IA\uj)^ and in turn (ogJ^{y) < +00. 

Theorems 3.3 and 3.4 and standard results of the theory of T-convergence, tell us that 31? low-energy 
sequences converge, up to subsequences, to a minimiser of the derived 2D model. This is the content of the 
following theorem. We refer the reader to (3.4) and the subsequent paragraph for the definition of the 3D 
total-energy functionals 

Theorem 3.5 (Splay-bend plate model). Suppose that the rescaled loads fh are such that fh/h? =" f weakly 
in L^(n,]R^) and satisfy the normalizing condition J^fi^dx = 0. Define the 2D total energy functional 
^SB ■ Aiso —R as 

■^ssiy) ■■= 4 T(y)- / f^^{x')-yix')dA, 

UJ 

where Agff is defined as in (3.21) and f^'^{x') := 4(^2 /(t', 2 : 3 ) dxa, for a.e. x' G uj. Suppose that {yh} is 
a low-energy sequence, viz. 

1 - ■^SBiVh) _ ,. infwi.2(n,R3) -^sb 


lim 
h^o h^ 


= lim 

h^o h^ 


=: m. 


Then, up to a subsequence, yh —> ysB in W^’^(f?,K^), where ysB G Aiso is a minimiser of the 2D model, 
that is 

^SBiVSB) = X^fsB- 

Moreover, m = ^sBivsB)- 

If we let / = 0 in the above theorem, we have 

min.^^ = min 4 ^ = 4 ™(ysB) = ^ (1 + 7)<5o^ 
and the associated fundamental form of ysB is given by ( 12 do/ 7 r^) diag(— 1 , 0 ). 
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Let US now fix a low-energy sequence {yh\ converging to a minimiser y G Aiso and rephrase the theorem in 
terms of the physical total energies defined in (2.6). Defining the deformations Vh{z',Z 3 ) = yh{z', z^/h) 
in the physical reference configuration we have ^gg{vfi)/h^ = min_ 4 ^^^ in view of (3.3). 

Equivalently, for a given small thickness ho, the approximate identity 


'^SbKo) - ^ y| Ay(a;') + ^diag(l,0) -k 7(Hy(x') + ^) jda;' 


+ fj,hl{l+ 7)^0 ( 

holds true, modulo terms of order higher than 3 in /iq. 


TT^ - 12 


yi-hlj f^^{x')-y{x')dx' (3.22) 


Remark 3.6. Clearly, the function Wq defined in (3.13) vanishes in SO(3). Also, by the standard inequality 
between arithmetic and geometric mean we have that jAp > 3(detF)^/^ for every F G with positive 

determinant, which proves that 

Wo(F) > y ipit) := t-l-logt, t>0. 

In particular, we have that Wo(^) = 0 iff F G SO(3) and that Wo{F) > C\F\'^ for every large |F|. Moreover, 
due to the regularity of Wo around SO(3), the energy density grows quadratically close to SO(3). These facts 
show that Wo satisfies Assumption 3.2. 


We now move to the twisted geometry. In this case, the (renormalized) spontaneous strain distribution is 
given by 

Ch,T{x3)=l+^ (^Nt{x3)<^Nt{x3)-^'^ |1R5^||oo =o(M, (3-23) 

where Nt is defined as in (3.5), and the (rescaled) stored energy density is 


W^ix3,F) 



F) ■ Cf^^j^ixs) - 3 - 2 log(det F) 


+ Wt,o;(det F) 


on every deformation gradient F G such that detF > 0. Proceeding similarly to the splay-bend case, 

we set Uh,T ■= \JCh,T, so that W^{x 3 ,F) = ITo(F17^^^(x 3)), being Wo defined as in (3.13). Note that, by 

Taylor-expanding W Ch,T around I, we get 


Uh^T =^ + h 





AIj' := Nx 


cos^/i isin( 2 /i) 0 \ 
isin( 2 /i) sinVi 0 j 
0 0 0 / 


where ||r^||oo = o{h), the positive constant 6 o is defined as in (3.15), and where 0:3 >->■ fiix^) given by (2.12) 
(with h = 1). Hence, we can write 


W^{x3,F) = Wo{F{l + hBl^ix3))), 


TDh _ 

.— 





(3.24) 


and we have that Bt in L°°((—1/2,1/2),, where Bt ■= —(5o (Mrixs) — |). Now, arguing as 

for the splay-bend case and using (3.17)-(3.18), we are left to derive (cfr (3.11)) the expression for Bt = 
—60 {Mt ~ where 


MtIxs) := {Nt{x3) 0 Nt(x3)) = 


cos^ (71/4-I-71x3/2) i sin (71/2-I-71x3) \ 

i sin (77/2-I-TTXs) sin^ (77/4-I-77x3/2) /’ 
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and to compute, for every G e 
1/2 


Q2(G') := J Q 2 {D + tG + BTit))dt 


- 1/2 


1/2 , 


= 2 fj, min 

DeSym(2) 


- 1/2 


D + t symG — Sq MT{t) + ^l 2 

O 


+ 7 tr^ ( D -\-tG — 6 q MT{t) + -^l 2 )> dt. 


The integrals 


/ (i + (i +1') = (i+ 


dt = —, 

TT 


- 1/2 


and 


give 


1/2 


t cos' 


(4 + 2 ') 


- 1/2 


1/2 


Mj' dt = 


-1/2 
1/2 

- 1/2 


1/2 I/tt 
I/tt 1/2 


-1/2 


1/2 

J t sin + Trt^ dt = 0, 


-1/2 


1/2 


— I/tt^ 0 


0 1/t, 


and J t Mt dt = 

-1/2 -1/2 

These computations, together with the fact that ti Mt = \Mt\ = 1, show that Q 2 t{G)/{2^) equals 


1 (isymGp + 7 tr^G) + ^ symG • diag(l, -1) + ^ 


where 


qr^D) := |Zlp + ytr^G - do 
It is easy to see that 


sym D 


1 2/tt 
2/7r 1 


--(l-7)trt^ 


min qxiD) = qx 

DeSym.(2) 


do/6 do/ir ] \ ^ 
do/TT do/6 \) V 


18 ^7r2j^0’ 


and in turn that 

Q2,t{^) = 2/1 


^(^IsymGp Tytr^G^ + ^symG • diag(l,-l) + (^ - ;^)^( 


Other straightforward computations show that, setting 


ar = 3 ^, 


At = i^diag(-l,l). 


^ — 47r^ — 48\ 

= M (--) 


>1 


(3.25) 


one has 

02,t(^) = ciT Q 2 [G — AtY + Pt, for every G G Sym(2). (3.26) 

To state the result pertaining to the twisted model, we define the functional • Aiso —^ [0,oo), where 
the class Aiso is defined in (3.8), as 


4 ™( 2 /) := 2 / Q2,T{Ay{x'))dx' 


JA 

12 


y I Ay(x')-i^diag(- 1 , 1 ) + 7Hg(a;')| dx'+ ^ (3-27) 
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We recall that ^ and 7 are the elastic constants appearing in (2.5) and defined in (3.18), respectively. As for 
the splay-bend case, well-known results of the theory of F-convergence easily imply the following theorem. 
We refer to (3.4) and to the subsequent paragraph for the notation related to the 3D models. 

Theorem 3.7 (Twisted plate model). Under the same assumptions on the family of (rescaled) loads {fh} as 
in Theorem 3.5, define : Aiso —> R as 


:= / r^{x')-y{x')Ax', 


(3.28) 


where 4*™ is given by (3.27) and /^™(a;') := 4 ^ 2 /( 2 ^^ 2 ^ 3 ) for a.e. x' S w. Suppose that {yu} is a 
low-energy sequence, viz. 

infyGwu2(n,B3) 

hm ——— = hm --i--^- =: m. 

h —^0 Jt h —^0 h 

Then, up to a subsequence, yt —> yr in W^’^(n,K^), where yx S Aiso is a minimiser of the 2D model, that 
is 


^•^^{yT) = min.^l*'". 

Aiso 

Moreover, m = 

In the case where the limiting load / is identically zero, we have that min^^^^ . 4 *"* = nainAiso 4*™ ^^id 
the minimisers of 4*™ 9-^6 given by the following lemma. 


Lemma 3.8. We have that 


mm 

A,so 


= 4-(yT) = ^ 

TT^ 


^lim _ jpHm 


12 


1 -|- 27 ^ ^ — 48 

1 -I- 7 


(3.29) 


where yx G Aiso is such that 

either Ay„. = diag ( — - 


12 ^n 


^(1 + 7 )’ 


or 


Ayj, = diag 0 , 


12So 


^(1 + 7) 


Proof. Clearly, a deformation y G Aiso which minimises the integrand of S’^"^{y) pointwise is a minimiser of 
£>hm |-] 2 g clagg Aiso- Seeking for such a minimiser and since det Ay(a;') = 0 a.e. in uj whenever y G Aiso, 
we consider the problem 

I |A — diag(a, —q;)|^-I- 7 tr^Aj = min^ |(1-|-7)(^-|-ry)^-I -2 a (a — ^-|- 


min 

AlGSym(2):det ^1=0 


where we have set a := —12ho/7i^ and used the notation 
Sym(2). Setting 

/(C,??) := (l + 7)(C + ’7)^ + 2a(a-C + ??) 

we have that the previous minimisation problem can be rewritten as min^^gR ming^(^ /(^, 77 ). Using the 
method of Lagrange multipliers it is then easy to check that 


e c 

C I 

and 


(3.30) 

to represent an arbitrary matrix A G 

5c(C.?7) := £.11-C, 


.’IC ) = >4) = fi£o >% ) = f 


2 f 1 + 27 

1 -I- 7 


where 


- b 


a 


2 Vl + 7 


^4 


1 -I- 7 


■4(2 


vf 


1 -I- 7 




a 


1 -f 7 


+ 4(2. 


Correspondingly, we have that the solutions to the minimum problem on the left hand side in (3.30) are 
A+ := diag (4,4) = diag 0^ and A" := diag (4, % ) = diag ^0, “ • 
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Now, since there exists y G Aiso such that Ky = A~^ or Ay = A (this corresponds to y{uj) being locally 
isometric to a cylinder), we have obtained that 


mm 

y&Aisa 


125o 


Ay{x') -^diag(-l,l) +^Ay{x')>dx' = 


cj mm 

AGSym(2):det A=0 


||A- diag(a,-a)!^ + 7 tr^A|, 


and in turn that 


pUm _ M I, ,1 fU+ ^+1 , - 47r2 - 48 


min4- = + 


y , 2/^1 +27 \ ~ 477^ — 48' 

hI“|“ (t+t'+^'—=— 


Substituting in the last expression the definition of a gives (3.29). 

□ 


Similarly to formula (3.22), the limiting plate theory for the twisted case can be expressed in terms of 
physical parameters by 


= 


fJ-hl 

12 


l2So 


+7Hy(a;')jdx' 

yh^SQ — 48 

^ 2 


1^1 - ^0 / /'™(a;') ■ y{x') dx', (3.31) 


for a given small thickness hp, where the approximate identity holds modulo terms of order higher than 3 in 
ho, and where y G Aiso and Vhg are a minimiser of the 2D model (3.27)-(3.28) and a low-energy (physical) 
deformation, respectively. 

To put into perspective the two plate models which we have derived for splay-bend and twisted nematic 
elastomer thin sheets, we conclude this section with a comparison with the case where a limiting plate model 
originates from a three-dimensional spontaneous strain distribution which is simpler, i.e., quadratic in the 
thickness variable. We see that, as expected, when the spontaneous strains are kinematically compatible, the 
limiting two-dimensional stored energy functional is minimised at the value zero. 

Remark 3.9 (The linear/quadratic case). Consider a system in the (physical) reference configuration ilfi 
endowed with a stored energy density of the form (2.5), with the spontaneous strain distribution given by 

Ch{z 3 ) = I + 6 o Z 3 P + rjo R, (3.32) 


for some constant and dimensionless symmetric matrices P and R. Moreover, Sq and ryo are real constants 
whose dimensions are inverse length and square of inverse length, respectively. Let us denote by Ch{x^) the 
rescaled spontaneous strain 


Chixs) := Chihxs) = I -I- ^0 hxs P + rjo {hx^f R, 

and by (x 3 ,F) >->• W^{x 3 ,F) the corresponding stored energy density. Defining Uh ■= we have 

W^{x 3 ,F) = Wo(i^l7^ ^(xa)) = lTo(F(I-|-hi?^(x 3 ))), withlTo defined as in (3.13) andR^(x 3 ) := (—i5o X 3 P/2-1- 
Th/h) Since B in L°°((—1/2,1/2),IR^^^), with B{x 3 ) := —S 0 X 3 P/ 2 , then Theorems 3.3 and 3.4 tell 

us that the limiting two-dimensional plate model is described by the energy functional 

J Q2 (yAyix') - ^P^ dx. 


, 1 

Aiso ^ y ^ 


(3.33) 














18 


V. AGOSTINIANI AND A. DESIMONE 


Indeed, for every G G we have 

1/2 

Qo(G) = 2a min 

DeSym(2) 

- 1/2 

1/2 

/• 

- 1/2 


D + t sym G — 


+ 7tr2 £, + iG- Uf 


= 2 a min 

DGSym(2) 


\Df +-ftr^D + t^ 


(5o - 2 


symG ——P 


■7tr^(' 




dt 


1/2 


'Q 2 {c 


t"Q 2 {G-"-^P]dt= ^Q.iG- 4 p]- 




- 1/2 


Note that the coefficient multiplying the purely quadratic term in (3.32) does not play any role. Referring 
to Subsection 2.2.3 for a discussion on the kinematic compatibility of (3.32), we observe that in each of the 
four cases where (3.32) is kinematically compatible, listed in the mentioned subsection, we have that P has at 
least one zero eigenvalue. Hence, we have that the functional (3.33) can be minimised to zero by a deformation 
y G Aiso such that y(w) is locally isometric to a plane when both the eigenvalues of P are zero, and such that 
y{oj) is locally a cylinder in all the other cases. 


4. Energy minimising shapes under zero loads 

The aim of this section is to give an explicit representation of the minimal energy configurations of the 
nematic sheets and to gain some physical insight on their behaviour. To do this, we start by characterising 
the deformations y realising the condition Ay = diag(fc, 0 ), for some constant fc 7 ^ 0 , under the constraint of 
being isometries. More explicitly, we look for a (smooth) deformation y : w —> such that 

{VyfV'y = I 2 , Ay = (^ J ^ , (4.1) 

or, equivalently, such that 

/ |9iyp diy ■ d2y \ _ ( I d \ f diy ■ div diy ■ d2v \ ( k 0 \ 

V i^iy • Say ISayp j 0 1 / ’ V Say • Siu ^ay • Sau / 0 0 / ’ 

where v := Siy A Say. It is easy to check that deformations y satisfying these conditions are defined up to 
arbitrary translations and superposed rotations. Hence, we will use the normalising conditions 

y(0,0)=0, Vy(0,0) = (ei|ea), (4.2) 

to construct one specific representative. 

Note that from the condition 9ay • div = 0 and from the identity (SiSay) • v + 9ay • div = 0, obtained 
by differentiating 9ay • u = 0 with respect to xi, one gets (SiSay) • u = 0. Moreover, by differentiating the 
conditions |Siyp = 1 and jSayp = 1 with respect to X 2 and cci, respectively, we obtain 

(SiSay) • Siy = (SiSay) • Say = 0. 

Hence, we have that Si Say = 0 in w. Similarly, using the condition 5ay ■ Sau = 0 and suitably differentiating 
the identities Say • u = 0, |Sayp = 1, and diy ■ Say = 0, one gets that SaSay = 0 in w. This fact, coupled with 
the information that the mixed derivatives of y vanish, says that y must be of the form 

y{xi,X2) = X2CPw{xi), 

for some c € and some smooth w : M —>• such that |c| = |w(a:i)| = 1 and c • w{xi) = 0 for every Xi, 

where we use the notation w = diw. Observe that V'y = {w \ c), that W = (z> | 0), with 

( C3W2 - C2W3 \ 

C1W3 - C3W1 
C 2 W 1 - C 1 W 2 / 
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and that condition diy ■ 82 ^ = 0 in now automatically satisfied. Note also that we have not exploited yet 
the information that diy ■ div = w ■ i> = k, which is going to determine the explicit expression of w. More 
precisely, the function w has to satisfy the following system of equations: 

< W • C = 0, 

\w ■ V = k. 


To proceed, we set f := w and choose c = 62 , so that the above system reduces to 

if? + fi = h 

\f 2 = 0 , 

I-/1/3 + /1/3 = k. 


Setting /i(a;i) = cos (d(a:i)) and fsixi) = sin (0(a;i)), we have that the first equation is satisfied, while the 
third equation reduces to 9{xi) = —fc, which yields 9{xi) = —kxi + 5, for some ^ G K. In the end, we have 
obtained that 

/ /i \ ( cos(-fca;i+^) \ / -^sm{-kxi + ^) + mi 

/= /2 = 0 ^ w = f = \ I cos(-A:a;i +0 + ^2 

\ fs J V sin(-fca;i +0 / \ m 3 

for some constant ^ G K, m G All in all, we have that if ?/ : a; —>■ is a smooth deformation satisfying 
(4.1), then d 2 y = c for some constant c G of unit length. Under the normalising assumption that c = 62 , 
the deformation y has the following expression 

/ -\s\a{-kxi + 0 \ 

y{xi,X 2 ) = X 2 & 2 + \ 0 +m, (4.3) 

\ I cos(-fca:i +0 / 

for some constants ^ G K and m G We can now choose ^ = 0 and m = (0, 0, —l/k)'^, so that 3iy(0, 0) = ei 
and y(0,0) = 0. Summarizing, the deformation 

T 

(4.4) 

fulfills condition (4.1) and the extra conditions (4.2). 

If we now look for some isometric deformation y realising the condition Ay = diag( 0 , k), for some constant 
k ^ 0, namely, such that 

iVyfVy = h, ^y=(o fc)> 

we can proceed similarly to the above and check that it must be of the form y(xi,X 2 ) = xic + w{x 2 ), for 
some c G and some smooth ic : K —)• such that |c| = \d 2 w{x 2 )\ = 1 and c • d 2 w{x 2 ) = 0 for every X 2 - 
Choosing c = ei, we easily arrive to the expression 

( 0 

y{xi,X2) = xiei+ \ lsm{kx2+0 
\ \ cos(fca;2 + C) 

Choosing ^ = 0 and m = (0,0, — 1/fc)^, we obtain the deformation 



y{xi,X2) = 


— sin(A;a;i), X 2 , y( cos(fcxi) — l) 
k k^ ' 


y{xi,X 2 ) 


Xi, i sin(A:a:2), ^(cos(fca;2) “ 


(4.6) 


fulfilling conditions (4.5) and the same extra conditions as y in (4.2). 

The spontaneous curvature exhibited by minimal energy configurations in twisted nematic elastomer sheets 
cannot be read off directly from the target curvature tensor. This is because the two-dimensional bending 
energy (3.27) cannot be minimised by minimising the integrand to zero, due to a geometric obstruction (there 
is no isometry of the plane with non-vanishing Gaussian curvature). This curvature is instead obtained by 
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solving a minimisation problem, as shown in Lemma 3.8. This lemma, coupled with the above discussion, 
says that the deformation y defined as in (4.4) with k = (a portion of a cylinder with axis parallel 

to the image through y of the line spanned by 62 , and with radius p := 7 r^(l + 7 )/( 12 ( 5 o)) and the deformation 
y defined as in (4.6) with k = (in this case, a portion of a cylinder with axis parallel to the image 

through y of the line spanned by ei, and with the same radius p) both realise the minimum for the 2D twisted 
energy functional. Nematic sheets with twisted texture are therefore bistable under zero loads, see Figure 2. 



Figure 2. Plot of minimal energy configurations for the free-energy functional defined 
in (3.27), arising from a twisted texture of the nematic director. Under zero loads, the system 
is stable in each of the two configurations, hence bistable. 


In the case of splay-bend textures, the curvature giving minimal energy can be predicted by simply reading 
it off from the target curvature tensor of the two-dimensional model (3.22) and therefore only deformations 
of type (4.4) with k = (a portion of a cylinder with axis parallel to the image through y of the line 

spanned by 62 , and with radius 7r^/(12^o)) are minimal energy states. This means that Gaussian curvature 
is suppressed in the splay-bend as well as in the twisted case, in the sense that the configurations exhibited 
by elastomer thin sheets in the absence of applied loads will be portions of cylindrical surfaces (with zero 
Gaussian curvature, as predicted in [35, 36] and observed experimentally in [30, 31, 34]). In both cases, these 
configurations carry non-zero residual stresses. In the twisted case, there will be also non-zero residual internal 
bending moments, due to the additional frustration caused by the non-attainability of the target curvature. In 
the splay-bend case the target curvature is attained, the bending energy is minimised to zero, and no residual 
moments arise. 

It is worth comparing the case of twist and splay-bend textures with a different scenario, in which the 
nematic director is kept constant along the thickness of the thin sheet, whereas the spontaneous strain ( 2 . 2 ) 
varies along the thickness through the magnitude parameter a. To the best of our knowledge, a system with 
these features has not yet been synthesized in a laboratory. At least in principle, this should be possible by 
realising a film with uniform alignment of the director perpendicular to the mid-suface (direction 63 ), and 
variable degree of order along the thickness {xs coordinate). 
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Using the notation of Subsection 2.1, let us suppose that the nematic director 'nP'{z^) is now constant, and 
equal to some n G S^, and that the (constant) parameter ah in (2.4) is now given by 

^h{zz) '■= 1 + ^3 € (~^/2,/i/2). 

ho 

The (physical) spontaneous strain of this system is therefore defined as 

Chizs) := d^^^(2:3) n (g) n + d)))^^^(z3)(l - n (g) n). 

Modelling the system using again the prototypical energy density (2.5), as in Subsection 3.1 we can define 
the rescaled energy densities Whix^,-), X 3 G (— 1 / 2 , 1 / 2 ), characterised by the (rescaled) spontaneous strains 

Chixs) := Chihxs) = [1 + j^hxsj n ® n + (^1 + 0:3j (l - n g) n) 

= l-2hB{x3) + R'^{x3), B{x3) := y ^ Q - n (g) , 

where ||i?^||oo = o{h). Proceeding as in Subsection 3.3, we obtain a limit 2D model whose free-energy 
functional is given by 


for every y G Aiso- Here, the 2x2 symmetric matrix M is given by the formula 


2 ho 


. - I 2 

(n g n) - y 


and (n g) n) is the 2x2 upper left part of n g n. In the case where n = 63, the spontaneous curvature tensor 
M reduces to 


M = 


Too 0 
0 mg 


mo := - 


OiO 

Qhn 


(4.8) 


Note that the Gaussian curvature associated with M is positive. However, as for the twisted case, where 
the spontaneous Gaussian curvature is negative, the observable minimal energy configurations will always 
exhibit zero Gaussian curvature (see Lemma 3.8), because of the isometry constraint they are subjected to. 
More precisely, some calculations show that every isometric deformation y G Aiso such that Ay = H+(s) or 
Ay = H_(s) for some s G [—fc/2, fc/2], with 


A±{s) := 




k^yj fc2—4 s2 
2 


k := mo 


/ 1 + 27 \ 

V 1 + 7 /’ 


(4.9) 


is such that 


yyj 12 


h 2 /1 + 27 




where the expression of 7 in terms of the 3D parameters is given in formula (3.18). Note that A_{—k/2) = 
H-|_(—fc/ 2 ) and A-(k/2) = A+(k/2), whereas H_(si) ^ H+(s 2 ) for every si,S 2 G (—fc/ 2 ,fc/ 2 ). Note also that 
the eigenvalues of H_|_(s) and H_(s) are always k and 0 for every s G [—fc/2,fc/2], and that 


H+(s) = i?+(s) 


k 0 
0 0 




H_(s) = R-{s) 


0 0 
0 k 




(4.10) 


where the columns of the rotation matrices R+{s) and R-{s) are the eigenvectors corresponding to k and 0 
and to 0 and k, respectively. The explicit expressions of R+{s) and R-{s) are the following: 


R±{s) = 


k + \/k'^ — 4s^ 
' 2 fc 


T 


± 


2s 


k+'\/ k^ — 4s^ 


k+^y k^—4s‘^ 
1 


(4.11) 
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In particular, for A^, the directions corresponding to the eigenvalue 0 are given, respectively, by the vector 
(1,1) in the case s = —kj2, by (0,1) in the case s = 0, and by (1,-1) in the case s = k/2. For As, the 
directions corresponding to the eigenvalue 0 are given, respectively, by the vector (1,1) in the case s = —k/2, 
by (1,0) in the case s = 0, and by (1, —1) in the case s = k/2. Therefore, through the matrices A+ and A- 
all the possible directions corresponding to the eigenvalue 0 (and in turn all the corresponding orthogonal 
eigenspaces corresponding to the eigenvector k) are represented. All in all, we have that 


s/ ■.= I A-|_(s) : s G 


k k 
2’ 2 


U < A_(s) : s G 


k k 
2’ 2 


|Adiag(A:,0)i?'^ : E G SO(2)| 


From the discussion leading to expression (4.3), we have that, given p > 0, the deformation defined as 


y(xi,X2) 



is an isometry such that Ay = diag(l/p, 0). Moreover, we have that 


y(0,0) = (0,0, p)^. 


Consider the rotation matrices 


:= 


sin a — cos a 0 \ 
cos a sin a 0 , 

0 0 1 / 


Vp(0,0) = (ei|e2). 


:= 


sin a — cos a 
cos a sin a 


(4.12) 


and note that Ea is a (7r/2 — Q;)-counterclockwise rotation taking the vector (cos a, sin a)^ into ( 0 , 1 )^. Now, 
setting p = 1/k, where k is defined as in ( 4 . 9 ), we define the deformation 


/ 


z/“(a;') := El o y o Ea{x') = 


sin ^ 


X\ Sin OL — X 2 cos Q. 

p 


i 


psm a sin 
cos a sin 

^ p cos a cos 

Simple computations show that € .4^50 and that, setting := A 92y", 


:.isina-:.2Cosa ^ ^ ^ gj^^ q,) 

+ sin a(xi cos a + X 2 sin a) 

) / 


xi Sin 0C—X2 cos a 
P 


Ay^ix') := (V'p“(:r'))^VV“(:r') = EI Ay{E^x') Ec. = ( 0 q ) ^ 


(4.13) 


Therefore, we have that Ayc = A_|_(s) for some s G [—fc/2, fc/2] and a i—>■ ?/“ is a continuous family of 
deformations minimising Note also that 1/“ satisfies the normalizing conditions (4.12), for every a. 

Figure 3 shows minimal energy deformed configurations obtained from the family a ^ y‘^. 

The existence of a continuous family of deformations with (constant) minimal energy shows that a nematic 
elastomer sheet with constant director (perpendicular to the mid-surface) and thickness-dependent magnitude 
of the spontaneous strain (this can be realised by varying the degree of the nematic order along the thickness) 
realises a “zero-stiffness” structure in the sense of [19]. These are structures that can undergo large elastic 
deformations without requiring external work. Figure 3 show that the nematic sheet can accomodate any 
level of twisting with negligible elastic energy in between two extreme states (a = 0 and a = 7r/2). Of course, 
zero-stiffness is an idealisation and, in a real system, effects that have not been taken into account in the 
model will lead to small, but non-zero loads in order to change shape. In the example of Sharon [23], edge 
effects cause energy storage which scales as h/^^'^. This is a higher scaling (with smaller stored energy in the 
thin film limit —>■ 0) with respect to the bending one {h^) that our dimensionally reduced theory is designed 
to resolve. As a consequence, the observed response is much “softer” than the one expected from the bending 
stiffness of a sheet. 

By contrast, sheets with twist texture are “bistable” in the sense of [20]: they exhibit two distinct possible 
stable shapes in the absence of loads (see Figure 2). Splay-bend sheets have only one shape minimising the 
energy under zero loads. 
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Figure 3. Plot of minimal energy configurations for the free-energy functional defined 
in (4.7)-(4.8), arising from constant director 63 along the thickness of the sheet. The con¬ 
figurations are elements of the continuous family of deformations (4.13). Under zero loads, 
the system has the same (minimal) energy in each of these configurations, hence it realises a 
structure with zero-stiffness to twisting. 
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